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Abstract 

We study domain growth in a nonlinear optical system useful to explore 
different scenarios that might occur in systems which do not relax to ther- 
modynamic equilibrium. Domains correspond to equivalent states of different 
circular polarization of light. We describe three dynamical regimes: a coars- 
ening regime in which dynamical scaling holds with a growth law dictated by 
curvature effects, a regime in which localized structures form, and a regime in 
which polarization domain walls are modulationally unstable and the system 
freezes in a labyrinthine pattern. 
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The problem of the growth of spatial domains of different phases has been thoroughly 
studied in the context of the dynamics of phase transitions: a system is placed in an un- 
stable state and one considers its relaxation to the state of thermodynamic equilibrium [|TJ . 
This process is dominated by the motion of domain walls and other defects. It is in this 
context that seminal ideas of selfsimilar evolution and dynamical scaling were introduced for 
nonequilibrium processes. Asymptotic domain growth laws, with their underlying physical 
mechanisms, have been well established, and dynamical scaling has been generally demon- 
strated. A growth law R(t) ~ t 1//2 holds for dynamics with no conservation law and domains 
made of equivalent phases. This law follows from the minimization of surface energy, and 
it has been shown to be robust against the appearance of point defects in systems with 
a discrete number of phases, three dimensional vortices or chiral domain walls 0. Other 
well known growth laws [Jj are R(t) ~ t 1//3 for systems with conserved order parameter and 
R(t) ~ t for nonconserved dynamics with a metastable phase ||, and also for hydrodynamic 
systems in spatial dimension d > 2 j4|. 

Domain growth in systems that do not approach a final state of thermodynamic equi- 
librium is much less understood. For example, the mechanisms underlying a growth law 
R(t) ~ t 1 / 5 in pattern forming systems in which the spatial coupling is non purely diffusive 
(Swift-Hohenberg equation) p have not been clearly identified. Other general issues that 
need to be considered are the role of hamiltonian vs. dissipative dynamics || , the effects 
of nonrelaxational dynamics such as one-dimensional motion of fronts between equivalent 
states and spiral formation 0, the emergence of localized structures (LS) |8|||, or transverse 



instabilities of domain walls leading to labyrinthine patterns [|10 . 

Driven nonlinear optical systems offer a wealth of opportunities for the study of pat- 
tern formation and other nonequilibrium processes in which the spatial coupling is caused 
by diffraction instead of diffusion. These systems are specially interesting because they 
naturally lead to the consideration of vectorial complex fields, being the vector character as- 
sociated with the polarization of light, and also because they often support the formation of 
LS [pTHl3[| . Such bright light spots are being actively considered for applications in parallel 



optical processing. Only very recently domain growth has been considered in some of these 
systems and some growth laws obtained from numerical simulations have been reported 
|T3|^r5||. However, clear mechanisms for the growth laws have often not been identified, and 
some of these laws do not correspond unambiguously to an asymptotic regime. In addition, 
the question of dynamical scaling has, in general, not been addressed. 

In this letter we consider a Kerr medium as a clear example of a nonlinear optical sys- 
tem in which many of the issues and scenarios mentioned above can be explored. We show 
that after switching-on a pump field, domain walls are formed which separate regions with 
different polarization of light. The dynamical evolution of these polarization domain walls 
leads to three different regimes. For high pump values there is a coarsening regime for which 
we demonstrate dynamical scaling with a growth law R(t) ~ t 1 / 2 . For lower pump values 
this process is contaminated by the emergence of LS formed by the collapse of polariza- 
tion domain walls to a stable bound structure. In a third regime the system evolves into 
a nearly frozen labyrinthine pattern caused by a transverse modulational instability of the 
polarization domain wall. These three qualitatively different regimes have been experimen- 
tally observed in another optical system [|HJ and considered in the realm of Swift-Hohenberg 
models NTH . 



Our calculations are based on a mean field model that describes the transverse spatio- 
temporal evolution of the two circularly polarized components of the electric field com- 
plex envelope, E + and E_, in an optical cavity filled with an isotropic self-defocusing Kerr 
medium and pumped with a linearly polarized real field E |17]18|| : 



d t E ± = -(1 - i9)E ± + iV\E ± + E -~i [\E±\ 2 + (3\E T \ 2 } E ± . (1) 

Here 9 is the cavity detuning, and is the laplacian in the transverse plane. Equations 
([[]) are damped and driven coupled Nonlinear Schrodinger equations which can be rewritten 

as 

dt E ± = -E ± -ij^, (2) 



where J-[E+, EJ\ is a real functional. Therefore, except for the linear dissipative term, 
the dynamics can be written in Hamiltonian form. This corresponds to a rather different 
dynamics than the normal relaxational dynamics considered in systems that approach a 
state of thermodynamic equilibrium. We will study different regimes for different values of 
the pump Eq. 

Equations (|l|) admit symmetric (I s+ = J s _) and asymmetric (I s+ ^ I s _) steady state 
homogeneous solutions, where I± = \E±\ 2 . The homogeneous symmetric solution is linearly 
stable for Eq < Eo a , while the asymmetric solutions only exist for Eq^ < Eqi < Eq and they 



are linearly stable for E 0j b < E 0:C < E |18|,|T9j . There are two equivalent homogeneous stable 
solutions for E 0>c < E , one in which I s+ 3> I s _ and the other one, obtained interchanging E + 
by E_, in which I s+ <C I s ^. These solutions are elliptically polarized, but very close to being 
circularly polarized, because one of the two circularly polarized components dominates. For 
simplicity we will call them the right and left circularly polarized solutions. If the pump field 
Eq is switched-on from Eq = to a value Eq > E ^ c , only the mode with zero wavenumber 
can initially grow from the initial condition E± = 0. One then expects that either of 
the two equivalent homogeneous solutions will locally grow and that domains separated 
by polarization walls will emerge. This is indeed the process that we study. We note, 
however, that a solution with a stripe pattern orthogonally polarized to the pump exists 



for Eq > E$ a [|T7|Jl8r| . This pattern solution is the one obtained by continuity from the 
homogeneous symmetric solution through a Turing-like instability. We have numerically 
checked that such solution remains stable for pump values Eq ^> E 0jC , but it is not the 
solution approached by the physical process just described of switching-on the pump to a 
value Eq > E 0tC . 

We find three different dynamical regimes for Eq > E Qc , summarized in fig. [I]. For 
Eq > Eq j2 domains grow and the system coarsens, for E 0i2 > E > _E 01 stable LS are 
formed, while for Eq^ > Eq > Eq^ a labyrinthine pattern emerges. These regimes are better 
understood by considering the evolution of an initial isolated polarization droplet: a circular 
domain of one of the solutions surrounded by the other solution. We find that the radius of 



the circular domain varies consistently with a curvature driven front motion. The normal 
front velocity v n (eikonal equation) follows a law of the form v n (r, t) = — / y(E )K(r, t), where 
k is the local curvature of the domain wall and j(Eq) is a coefficient that depends on the 
pump field amplitude. For a circular domain we get dR(t)/dt = —^(Eo)/R{t). In figure [I] 
we show the function 7(i?o) as obtained from the numerical solution of eqs. in a two- 
dimensional system for relatively large initial droplets. Notice that 7(-E'o) changes sign at 
Eq = Eo t i, which indicates a change from droplet shrinkage to droplet growth. 

We first consider the regime of domain coarsening which occurs for Eq > Eq^- In this 
regime 7(-Eo) > and an isolated drop shrinks to zero radius. In the general dynamics start- 
ing from random initial conditions around E± = 0, sharp domain walls are initially formed 
and they evolve reducing their curvature. The system approaches a final homogeneous state 
in which one of the two circularly polarized solutions fills the whole system. In order to 
characterize the coarsening process we have calculated the pair correlation function of I + 
and defined as C/ ± (r, t) = (i±(x + r, i)J ± (x, £)}. The average (. . .) is performed over the 
set of points x (and additionally over a set of 100 different random initial conditions). Due 
to the symmetry of the problem Cj + = Cj_ = C. Results for the circularly averaged corre- 
lation function C(r,t) are shown in fig. ^. The mean size L(t) of the domains is calculated 
as the distance at which C(r,t) takes half its value at the origin, i.e., C(L(t),t) = |C(0, t). 
We obtain a well defined asymptotic growth law Lit) ~ t 1 ! 2 that follows from domain wall 
motion driven by curvature effects. We have further obtained that the dynamics is self- 
similar, i.e., that there is dynamical scaling. This is seen in fig. [| where we plot C{r,t) 
before and after rescaling the spatial coordinate of the system with the characteristic do- 
main size L{t). We observe that curves for different times in the scaling regime collapse to 
the single scaling function after rescaling. These results coincide with those obtained for 
many thermodynamic systems with nonconserved order parameter We note, however, 
that in our case the dynamics does not follow the minimization of any obvious energy and 
that surface tension is not a proper concept for the diffractive spatial coupling considered 
in optical systems. 
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We next address the regime of formation of LS (i?o,2 > Eq > Eo,i). In this regime, as 
in the previous case, 7(-Eq) > 0, and a large isolated droplet initially shrinks with a radius 
decreasing as R(t) ~ i -1 / 2 . However the shrinkage stops at a well defined final value of 
the radius. Initial droplets with a smaller radius grow to this final stable radius. In the 
general dynamics following the switch-on of the pump, domain walls are initially formed. 
They first evolve reducing their length as in the coarsening regime. But while in that regime 
a closed loop disappears, here it collapses to a stable LS formed by a bound state of the 
domain wall. The final state is composed of stretched domain walls and LS. To understand 
this process is convenient to consider the form of the polarization domain walls in a d — 1 
geometry, as shown in figure [3| An isolated d = 1 domain wall is stationary. We observe 
that the intensity profiles of the walls do not approach monotonically the asymptotic value 
of the homogeneous state. When several domain walls are created in the transient dynamics, 
they interact with each other. Since the front profiles have oscillatory tails, the interaction 
between two walls can lead to repulsive forces ||. As a consequence, LS formed by bound 
domain walls can be formed which stop the coarsening process. These oscillatory tails are 
less important the larger is E (see figure |3|). However, we find that for all the values of 
E which we have explored (up to E = 10), this effect is enough to stop coarsening in 
d = 1: a frozen pattern state is always dynamically reached p0| . What happens in our 



d = 2 situation is a competition between the d = 1 repulsive effect between walls and the 
curvature effect that tends to reduce a droplet to zero radius. When the repulsive force is 
large enough, it might counterbalance the shrinkage process driven by curvature, and thus 
leads to the formation of a LS. This happens for E ^ < E Q < E ^- The mechanism is the one 
also discussed in [fH5| . These structures can be seen as a hole of I + (J_) in the background 
of a circularly + - polarized ( — polarized) state, together with a peak of 7_ {!+)■ Since the 
oscillatory tails are larger as E decreases, the size of the LS decreases with E . We have 
found a perfect linear dependence of the radius of the LS with Eq. In figure |] we show a 
plot of a LS together with its transverse profile. Note that the intensity in the LS is greater 
than in the surrounding background. 
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We finally discuss the regime of labyrinthine pattern formation which occurs for Eq < 
E 0j i : switching-on the pump produces a very dense pattern of domain walls that repel each 
other. In this regime j(Eq) < 0, and an isolated droplet of arbitrary small size grows as 
R(t) ~ t 1 ! 2 . In an infinitely large system the droplet would grow without limit, but with 
periodic boundary conditions it grows until the domain wall interacts with itself. Repulsion 
of the domain wall leads to a labyrinthine pattern as shown in fig. [5]. An independent way 
of identifying the value Eq = £7o,i> below which labyrinthine patterns emerge, is by a linear 
stability analysis in d = 2 of the d = 1 domain wall profile. We have numerically obtained 
that such flat domain wall has a transverse modulational instability for values of the pump 
amplitude for which ^(Eq) < 0. We find a longwavelength instability in which arbitrary 
small wavenumbers become unstable for E < E i (see fig. ^). This is reminiscent of the 
situation described for vectorial Second Harmonic Generation fl5|| . In physical terms, both 
the droplet growth and the modulational instability indicate that the system prefers to have 
the longest possible domain walls, or equivalently the largest possible curvature. This leads 
to a nearly frozen state in which the oscillatory tails of the domain walls prevent their self- 
crossing and in which coarsening is suppressed. LS might form, but their natural tendency 
to grow is stopped by surrounding walls. 

In summary, we have described a situation in nonlinear optics in which many of the 
generic issues and possible scenarios of domain growth in nont her mo dynamic systems occur. 
In spite of the nonrelaxational dynamics we have found a regime of selfsimilar evolution with 
a growth law characteristic of curvature driven motion. In other regimes, obtained just by 
changing the pump amplitude, domain growth is contaminated by the emergence of LS or 
suppressed by an instability of the domain wall that leads to a nearly frozen labyrinthine 
pattern. Domain walls and LS are here associated with the polarization vectorial degree of 
freedom of light. 
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FIGURES 
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FIG. 1. Coefficient 7(^0) as denned in the text. Snapshots of typical configurations for the 
total intensity I = i+ + I- at late times are shown for each of the three dynamical regimes. The 
vertical lines identify the values of £0,1 = 1.552 and £0,2 = 1-700. 
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FIG. 2. a) Spherical averaged correlation function for Eq = 1.8 > £0,2 a * times t = 64.8, 91.8, 
11.8, 145.8 and 172.8 and b) after scaling r with the domain size Lj + {t). The inset shows the 
domain growth law LiAt) ~ t 1 ! 2 . 






FIG. 3. 
E = 1.2 < 



10 20 30 40 10 20 30 40 

d = 1 intensity profiles of the polarization domain wall for two pump values: 
^0,1) Eq = 2> Eq$- 
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FIG. 4. Total field intensity (I + + I ) of a LS and transverse profile of I + and i_ for = 1-6- 
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FIG. 6. Transverse modulational instability for a flat domain wall. E = 1.4. Snapshots at 
times t = 0, 400, 600, 860. 
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